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(54) Title: IMAGING METHOD AND DEVICE USING SHEARING WAVES 

(54) Titre: PROCEDE ET DISPOSITIF D'IMAGERIE UTILISANT LES ONDES DE CISAILLEMENT 
(57) Abstract 

The invention concerns an imaging method 
for observing the propagation of a low-frequency 
shearing pulse wave simultaneously in multiple 
points of a diffusing viscoelastic medium (1). The 
method consists in transmitting at very high rate 
ultrasonic compression waves enabling to obtain a 
succession of images of the medium; then in delayed 
processing of the resulting images by intercorrelation 
to determine in each point of each image the 
movements of the medium while the shearing wave 
Is being propagated. 

(57) Abrege- 

Proc&te d'imagerie pour oberver la propaga- 
tion d'une onde impulsionnelle de cisaillement basse 
frequence simul landmen t en une multitude de points 
d'un milieu viscoelastique diffusant (I). A cet effet, 
on 6met a cadence ultra rapide des ondes ultrasonores 
de compression qui permettent d'obtenir une suc- 
cession d* images du milieu, puis on traite en temps 
Clifford les images ainsi obtenues par intercorrelation, 

pour determiner en chaque point de chaque image osp _ qjqual signal processor 

les mouvements du milieu lors de la propagaUon de cpu « central processing signal 

I 'onde de cisaillement. <* - GgiEjura 

S Z ADDER CIRCUIT 



Imaging method and device using shear waves. 

The present invention relates to imaging methods and 
devices using shear waves. 

5 

More particularly, the invention relates to a method of 
imaging using shear waves to observe a diffusing 
viscoelastic medium which contains particles reflecting 
ultrasonic compression waves, in which method an 
10 elastic shear wave is generated in the viscoelastic 
medium and the displacement of the viscoelastic medium 
subjected to said shear wave is observed by means of at 
least one ultrasonic compression wave. 

15 Document US-A-5 810 731 describes an example of such a 
method, in which the shear wave is generated locally 
inside the observed viscoelastic medium, by means of 
the radiation pressure of a modulated ultrasound wave 
focussed on a point to be observed. An additional 

20 ultrasound wave is then dispatched to this focal point, 
the reflection of the wave making it possible to 
ascertain certain propagation parameters of the shear 
wave (in particular the dynamic viscosity of the medium 
and its shear modulus) in the vicinity of the 

25 abovementioned focal point. 

This technique has the drawback of allowing the 
analysis of just a single point of the viscoelastic 
medium under study each time a shear wave is generated. 
3 0 If one wishes to obtain a complete image of the 
observed viscoelastic medium, it is necessary to repeat 
the operation a very large number of times, this 
involving a considerable idle time (for example, 
several minutes) to obtain this image. 

35 

This considerable idle time renders this prior art 
method impractical to use. 

Moreover, such an idle time may impede the use of said 


2 

method to obtain an image of a living tissue, which is always in motion. 


The above discussion of background art is included to explain the 
context of the present invention. It is not to be taken as an admission that any 
5 of the documents or other material referred to was published, known or part of 
the common general knowledge in Australia at the priority date of any one of 
the claims of this specification. 

An aim of the present invention is in particular to alleviate these 
10 drawbacks. 

Viewed from one aspect, the present invention provides a method of 
imaging using shear waves to observe a diffusing viscoelastic medium which 
contains particles reflecting ultrasonic compression waves, in which method 

15 an elastic shear wave is generated in the viscoelastic medium and the 
displacement of the viscoelastic medium subjected to said shear wave is 
observed by means of at least one ultrasonic compression wave, 

wherein the shear wave is generated by applying to the viscoelastic 
medium an excitation having the form of a low-frequency pulse which 

20 exhibits a central frequency f of between 20 and 5000 Hz, this low- frequency 
pulse exhibiting a duration of between l/2f and 20/f, 

wherein it comprises a propagation observation step in the course of which the 
propagation of the shear wave is observed simultaneously at a multitude of 
points in the observed medium, these points forming a substantially 
25 continuous observation field extending at least along a first axis, this shear 
wave propagation observation step consisting in: 

emitting into the observed medium a succession of at least 10 shots of 
ultrasonic compression waves at a rate of between 100 and 100 000 shots per 
second, 

30 - detecting and recording in real time the echoes generated by the 
reflecting particles of the viscoelastic medium at each ultrasonic wave shot, 


2a 

these echoes corresponding (directly or indirectly) to successive images of the 
observed medium, and 

wherein said method furthermore comprises a subsequent image processing 
step in the course of which the images thus obtained are processed at a later 
time at least by cross-correlation between successive images, so as to 
determine at each point of the observation field a motion parameter chosen 
between the displacement and the strain of the viscoelastic medium, in such a 
way as to thus obtain a succession of images showing the evolution of the 
motion parameter of the viscoelastic medium under the effect of the 
propagation 
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of the shear wave. 

By virtue of these arrangements, a film is obtained 
which clearly illustrates the propagation of the shear 
5 wave in the viscoelastic medium, which may make it 
possible for example, in medical applications, to 
directly tag cancerous areas in the tissues of a 
patient: the propagation of the shear waves in fact 
occurs there very differently from the neighboring 
10 areas. 

This tagging is performed much more easily than through 
conventional observation by simple ultrasound 
echography, since the propagation of the shear waves is 

15 dependent on the shear modulus of the medium, itself 
varying greatly between an area of healthy tissues and 
an area of cancerous tissues : the shear modulus 
typically varies in a ratio of 1 to 30 between a 
healthy area and a cancerous area, whereas the bulk 

20 modulus, which governs the propagation of the acoustic 
compression waves used in ultrasound echography, varies 
by only the order of 5% between a healthy tissue and a 
cancerous tissue. 

25 It will be noted that the film obtained illustrates the 
propagation of the shear wave much more clearly than 
the simple succession of images given by the reflecting 
particles of the medium, since said film makes it 
possible to view at each instant the areas of the 

30 observed medium which undergo motions of the same 
magnitude on account of the propagation of the shear 
wave, whereas the succession of images of the 
reflecting particles would make it possible to view 
only a haze of bright points in motion. 

35 

In preferred embodiments of the method according to the 
invention, recourse may furthermore possibly be had to 
one and/or to the other of the following arrangements: 
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the duration of the low- frequency pulse is 
between l/2f and 2f; 

the central frequency of the low- frequency 
pulse is between 30 and 1000 Hz; 
5 - the observed viscoelastic medium consists of a 

living body comprising at least one internal organ 
subjected to pulsatile motions, the low- frequency pulse 
which generates the shear wave being constituted by a 
pulsatile motion of said internal organ; 
10 - the observed viscoelastic medium is delimited 

by an outside surface and the low- frequency pulse is 
applied in the vicinity of this outside surface; 

the low- frequency pulse is applied by a means 
of excitation chosen from: 
15 • an acoustic wave generated by at least one 

acoustic transducer, 

• and a shock generated locally by physical 

contact in the vicinity of the outside surface of the 
viscoelastic medium; 

2 0 - the ultrasonic compression wave shots are 

emitted and the echoes generated by the reflecting 
particles of the viscoelastic medium are detected by 
means of a bank of transducers which comprises at least 
one transducer and which is arranged in contact with 
25 the outside surface of the viscoelastic medium, the 
shear wave being applied to the viscoelastic medium by 
imposing a pulsatile displacement on said bank of 
transducers ,- 

said motion parameter is the strain of the 
30 viscoelastic medium: this arrangement is particularly 
useful in the last case envisaged hereinabove, since it 
makes it possible to dispense with the displacement of 
the bank of transducers, which displacement would 
otherwise disturb the measurement of the displacement 

3 5 of the points of the observation field; 

in the course of the observation of the 
propagation of the shear wave, between 100 and 10 000 
ultrasonic compression wave shots are emitted at a rate 
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of between 100 and 100 000 shots per second; 

the observation field extends at least along a 
plane comprising on the one hand, the first axis and on 
the other hand, a second axis perpendicular to the 
5 first axis; 

in the course of the propagation observation 
step, a bank of several acoustic transducers arranged 
at least along the second axis is used to emit the 
ultrasonic compression wave shots and detect the echoes 

10 generated by the reflecting particles of the 
viscoelastic medium, the echoes detected by each 
acoustic transducer being stored directly without prior 
processing in the course of the propagation observation 
step, and the image processing step comprising a 

15 preliminary substep of forming paths in the course of 
which an image of the viscoelastic medium corresponding 
to each ultrasonic compression wave shot is generated 
by combining at least some of the echoes received by 
the various transducers; 

2 0 - the image processing step is followed 

(immediately or otherwise) by a viewing step in the 
course of which a film consisting of the succession of 
processed images is viewed under slow motion, each 
point of each image exhibiting an optical parameter 

25 which varies according to the value of the motion 
parameter assigned to this point; 

the optical parameter is chosen from the gray 
level and the chromatic level; 

the image processing step is followed 

30 (immediately or otherwise) by a mapping step in the 
course of which, on the basis of the evolution of the 
motion parameter over the course of time in the 
observation field, at least one propagation parameter 
is calculated for the shear wave at at least some 

35 points of the observation field; 

the shear wave propagation parameter which is 
calculated in the course of the mapping step is chosen 
from the speed of the shear waves, the shear modulus, 


6 

the attenuation of the shear waves, the shear elasticity, and the shear viscosity. 


Viewed from another aspect, the present invention provides an imaging 
device using shear waves to observe a diffusing viscoelastic medium which 
5 contains particles reflecting ultrasonic compression waves, this device 
comprising means of excitation for generating an elastic shear wave in the 
viscoelastic medium and means of acquisition for observing, by means of at 
least one ultrasonic compression wave, the displacement of the viscoelastic 
medium subjected to said shear wave, 
„ „ 10 wherein the means of excitation are adapted to apply to the viscoelastic 

: medium an excitation having the form of a low- frequency pulse which 

• • • 

exhibits a central frequency f of between 20 and 5000 Hz, this low-frequency 

• • • 

j pulse exhibiting a duration of between l/2f and 20/f, and 

wherein the means of acquisition are adapted to observe the propagation of the 
15 shear wave simultaneously at a multitude of points in the observed medium, 

„ . these points forming a substantially continuous observation field extending at 

• • • 

ry\ least along a first axis, said means of acquisition being adapted for: 

- emitting into the observed medium a succession of at least 10 

• • • 

shots of ultrasonic compression waves at a rate of between 100 and 100 000 
20 shots per second, 

• • 

: * • - detecting and recording in real time the echoes generated by the 

reflecting particles of the viscoelastic medium at each ultrasonic wave shot, 
these echoes corresponding (directly or indirectly) to successive images of the 
observed medium, and 
25 wherein said device furthermore comprises image processing means adapted 
to process at a later time the images obtained by the observation means, at 
least by cross-correlation between successive images, so as to determine at 
each point of the observation field a motion parameter chosen between the 
displacement and the strain of the viscoelastic medium, in such a way as 


30 
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to thus obtain a succession of images showing the 
evolution of the motion parameter of the viscoelastic 
medium under the effect of the propagation of the shear 
wave . 

5 

In preferred embodiments of the device according to the 
invention, recourse may moreover possibly be had to one 
and/or to the other of the following arrangements: 

the means of observation comprise a bank of 

10 transducers which includes at least one transducer and 
which is adapted to be arranged in contact with an 
outside surface delimiting the viscoelastic medium, the 
means of excitation being adapted to impose a pulsatile 
displacement on said bank of transducers; 

15 - said motion parameter is the strain of the 

viscoelastic medium; 

- the observation field extends at least along a 
plane comprising on the one hand, the first axis and on 
the other hand, a second axis perpendicular to the 

20 first axis, the bank of transducers comprising several 
transducers arranged at least along the second axis, 
control means being provided for selectively operating 
the device either in the mode of imaging by shear 
waves, or in a standard echography mode making it 

25 possible to acquire between 10 and 100 images per 
second. 

Other characteristics and advantages of the invention 
will become apparent in the course of the following 
30 description of several of its embodiments, given by way 
of nonlimiting examples, in conjunction with the 
appended drawings. 

In the drawings: 
35 - figure 1 is a schematic view of a device for 

imaging by shear waves according to one embodiment of 
the invention, 

figure 2 is a detail view showing a variant of 


WO 00/55616 


- 8 - 


PO/FROO/00599 


the device of figure 1. 

In the various figures, the same references designate 
identical or similar elements. 

5 

Figure 1 represents an example of a device for imaging 
by shear waves according to the invention, for studying 
the propagation of elastic shear waves in a 
viscoelastic medium 1 which is diffusing with regard to 
10 ultrasonic compression waves, and which may for example 
be : 

an inert body, in particular in the case of 
quality control for industrial applications, in 
particular agro-food applications, 
15 or a living body, for example a part of the 

body of a patient, in the case of medical applications. 

This device comprises an acoustic transducer such as a 
loudspeaker 2 or a vibrator can which is arranged 
20 against the outside surface 3 of the observed medium 1, 
this surface 3 consisting for example of the skin of 
the patient in medical applications. 

The loudspeaker 2 can be controlled by a microcomputer 
25 4, for example by way of a low-frequency pulse 
generator circuit G (this circuit can consist in 
particular of the sound card of the microcomputer 4) 
and of an amplifier A, so as to apply an excitation in 
the form of a low- frequency pulse to the surface 3 of 
3 0 the observed medium, so as to generate a shear wave in 
the viscoelastic medium 1. 

This low- frequency pulse generally exhibits an 
amplitude which may be of the order of 1 mm and a 
35 central frequency f of between 2 0 and 5000 Hz, applied 
for a duration of between 1/2 f and 20/f. Preferably, 
the duration of application of the low-frequency pulse 
is between l/2f and 2f and the frequency f is between 
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30 and 1000 Hz, this frequency typically being of the 
order of 50 Hz. 

As a variant, the acoustic shear wave could also be 
5 obtained (with the abovementioned amplitude and 
frequency characteristics) : 

- via a shock generated locally by physical 
contact in the vicinity of the outside surface of the 
viscoelastic medium, via at least one automatic 
10 mechanical actuator controlled by the microcomputer 4, 

via a shock applied manually by physical 
contact in the vicinity of the outside surface of the 
viscoelastic medium, 

or else, via a natural pulsatile motion of an 
15 internal organ of the human or animal body (for example 
a beating of the heart) in medical applications. 

The elastic shear wave produced by the loudspeaker 2 
moves with a relatively low speed Cs, of the order of a 
20 few m/s (typically, from 1 to 10 m/s in the human 
body) , producing internal motions in the observed 
viscoelastic medium 1. 

These motions are followed by dispatching into the 
25 medium 1, ultrasonic compression waves which interact 
with the diffusing particles 5 contained in the medium 
1, which particles are reflecting in respect of 
ultrasonic compression waves. The particles 5 may be 
constituted by any heterogeneity of the medium 1, and 
30 in particular, when dealing with a medical application, 
by collagen particles present in human tissues. 

To observe the propagation of the shear wave, use is 
therefore made of an ultrasound probe 6 arranged 
35 against the outside surface 3 of the observed medium 1. 
This probe dispatches, along an axis X, pulses of 
ultrasonic compression waves of the type of those 
commonly used in echography, at a frequency of for 


WO 00/55616 


- 10 - 


PO/FR00/00599 


example between 1 and 100 MHz and preferably between 3 
and 15 MHz. It will be noted that the probe 6 can be 
arranged : 

- either on the same side of the medium 1 as the 
5 loudspeaker 2, as represented in figure 1, 

- or on the opposite side from the loudspeaker 2 
with respect to the medium 1, 

- or in any other position, for example in a 
transverse arrangement with respect to the loudspeaker 

10 2. 

The ultrasound probe 6 consists of a bank of n 
ultrasound transducers Tl, T2 , . . . , Ti, . . . , Tn, n 
being an integer number n at least equal to 1. 

15 

This probe 6 usually takes the form of a linear strip 
which can comprise for example n = 128 transducers 
aligned along an axis Y perpendicular to the axis X, 
simultaneously dispatching their ultrasound wave pulses 
20 in such a way as to generate a "plane" wave (that is to 
say in this instance a wave whose wave front is 
straight in the X, Y plane) or any other type of wave 
illuminating the entire observation field. 

25 As a variant, the bank 2 of transducers may possibly be 
reduced to a single transducer Tl, or conversely take 
the form of a two-dimensional array extending for 
example along a plane perpendicular to the X axis. 

30 Each of the transducers Tl, T2, ... Tn is controlled by 
the microcomputer 4 or by a central processing unit CPU 
(which is contained for example in an electronic rack 7 
linked by a flexible cable to the probe 6) , so as to 
emit successive shots of ultrasonic compression waves 

35 into the medium 2, in the course of an observation 
phase which may last for example less than one second 
and in the course of which p shots of ultrasonic 
compression waves are emitted (p being an integer lying 
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between 100 and 10 000 and preferably between 1000 and 
100 000), at a rate of between 100 and 100 000 shots 
per second and preferably between 1000 and 100 000 
shots per second, in particular between 1000 and 10 000 
5 shots per second (this rate is limited by the time for 
the outward/return journey of the compression wave in 
the medium 1, hence by the thickness of the medium l in 
the direction X: it is in fact necessary for all the 
echoes generated by the compression wave to have been 
10 received by the probe 6 before a new compression wave 
is dispatched) . 

The ultrasonic compression wave shots of the 
observation phase preferably begin just before the 
15 emission of the shear wave. 

Moreover, in the case where the shear wave is generated 
by a pulsatile motion of an organ of a living body, it 
is advantageously possible to synchronize the starting 

2 0 of the ultrasonic compression wave shots with this 

pulsatile motion. For example in the case where the 
shear wave is generated by a heartbeat, it is possible 
to synchronize the starting of the ultrasonic 
compression wave shots with a chosen phase of the 
25 electrocardiogram. 

Each of these shots gives rise to the propagation of an 
ultrasonic compression wave in the medium 1, with a 
much higher speed of propagation than the shear waves, 

3 0 for example of the order of 1500 m/s in the human body. 

The ultrasonic wave thus generated interacts with the 
reflecting particles 5, thereby generating echoes or 
other similar disturbances of the signal, which are 
35 known per se as "speckle noise" within the field of 
echography. 

This "speckle noise" is picked up by the transducers 
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Tl, Tn after each shot. The signal sij (t) thus 

picked up by each transducer Ti after shot No. j is 
firstly sampled at high frequency (for example from 3 0 
to 100 MHz) and digitized in real time (for example on 
5 8 bits, or in certain cases on 1 bit) by a sampler 
belonging to the rack 7 and linked to this transducer, 
respectively El, E2, ... En. 

The signal sij (t) thus sampled and digitized is then 
10 stored, also in real time, in a memory Mi belonging to 
the rack 7 and specific to the transducer Ti . 

Each memory Mi exhibits for example a capacity of the 
order of 1 Mb, and contains all the signals sij (t) 
15 received successively for shots j = 1 to p. 

At a later time, after the storage of all the signals 
sij (t) corresponding to one and the same shear wave 
propagation, the central processing unit CPU has these 
20 signals reprocessed by a summator circuit S belonging 
to the rack 7 (or else it performs this processing 
itself, or else said processing can be performed in the 
microcomputer 4) , through a conventional process of 
path formation. 

25 

Signals Sj (x, y) are thus generated, each 
corresponding to the image of the observation field 
after shot No. j . 

3 0 For example, a signal Sj (t) can be determined through 
the following formula: 

n 

Sj(t) = £ a, (x, y).sij[t(x, y) + d, (x, y) / V] 
i-i 

35 where: 

sij is the raw signal sensed by transducer No. 
i after ultrasonic compression wave shot No. j, 
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t(x,y) is the time taken by the ultrasonic 
compression wave to reach the point of the observation 
field with coordinates (x,y), with t = 0 at the start 
of shot No. j , 

5 - di(x,y) is the distance between the point of 

the observation field with coordinates (x,y) and 
transducer No. i, or an approximation of this distance, 

V is the average speed of propagation of the 
ultrasonic acoustic compression waves in the observed 
10 viscoelastic medium, 

and oti(x,y) is a weighting coefficient allowing 
for apodization laws (in practice, in many cases it 
will be possible to consider that cti(x,y)=l) . 

15 The above formula applies mutatis mutandis when the 
observation field is 3 -dimensional (plane array of 
transducers) , on replacing the spatial coordinates 
(x,y) by (x,y, z) . 

20 When the probe 6 comprises just one transducer, the 
path formation step is unnecessary, and we have 
directly Sj (x) = sj[2.x/V], with the same notation as 
above . 

25 After the path formation step, should there be one, the 
central processing unit CPU stores in a central memory 
M belonging to the rack 7 the image signals Sj(x,y) or 
Sj (x) or Sj(x,y, z), which each correspond to shot 
No. j. These signals may also be stored in the 

30 microcomputer 4 when it performs the image processing 
itself. 

These images are then processed in pairs, again at a 
later time, by cross-correlation. This cross- 
35 correlation can be carried out in a DSP circuit 
belonging to the rack 7, or be programmed into the 
central processing unit CPU or into the microcomputer 
4 . 
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By way of example, this cross-correlation can be done 
by comparing the signals Sj (x,y) and Sj + 1 (x,y) (in 
the case of a two-dimensional observation field) over 
5 sliding spatial windows of predetermined length Ax, 
which may range for example from X to 10X, where X is 
the wavelength of the ultrasonic compression waves 
(i.e. around 0.42 to 4.2 mm at 3.5 MHz in water or the 
human body) . Moreover, the abovementioned windows may 
10 overlap one another over around 2 0% of their length 
along the X axis. 

In the course of this cross-correlation process, a 
cross-correlation function <Sj(x,y), Sj+l(x,y)> is 
15 maximized so as to determine the displacement undergone 
by each particle 5 giving rise to an ultrasound echo, 
in the direction X. 

Examples of such cross-correlation calculations are 
20 given in the state of the art, in particular by 
O'Donnell et al . ("Internal displacement and strain 
imaging using speckle tracking", IEEE transactions on 
ultrasonic, f erroelectrics, and frequency control, 
vol. 41, No. 3, May 1994, p. 314-325) and by Ophir et 
25 al. ("Elastography : a quantitative method for imaging 
the elasticity of biological tissues", Ultrasonic 
imag., vol. 13, p. 111-134, 1991). 

A succession of displacement fields Djx(x,y) of the 
30 medium 1 under the effect of the shear wave, in the 
direction X, is thus obtained. 

This succession of displacement fields is stored in the 
memory M or in the microcomputer 4 and can be viewed, 
35 in particular, by means of the screen 4a of the 
microcomputer, in the form of a slow-motion film where 
the value of the displacements is illustrated by an 
optical parameter such as by a gray level or by a 
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chromatic level. 

The differences in propagation of the shear wave 
between the areas of different characteristics of the 
5 medium 1, for example healthy tissues and cancerous 
tissues in the case of a medical application, are thus 
viewed perfectly. 

This shear wave propagation film can moreover be 
10 superimposed with a conventional echography image, 
which may be generated by the device described above, 
which is able to operate: 

either in the shear wave imaging mode, 
or in a standard echography mode, as a function 
15 of the commands received for example from the keyboard 
4b of the microcomputer. 

Moreover, it is also possible to calculate, rather than 
the displacements of each point of the observed medium 
20 1, the strains Ejx(x,y) of the medium in the direction 
X, that is to say the derivatives of the displacements 
Djx(x,y) with respect to x. 

These successive strain fields are usable as before to 
25 clearly view the propagation of the shear wave in the 
form of a film, and furthermore exhibit the advantage 
of dispensing with the displacements of the probe 6 
with respect to the observed medium 1 . 

30 This variant is especially beneficial in the embodiment 
of figure 2, where the probe 6 is carried by the 
loudspeaker or vibrator 2, this necessarily involving 
motions of said probe since it is the latter which 
itself then generates the shear wave. 

35 

On the basis of the displacement or strain fields, it 
is possible, if appropriate, to proceed moreover to a 
mapping step in the course of which, on the basis of 
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the evolution of the motion parameter (displacement or 
strain) over time in the observation field X, Y (or X 
in the case of a single transducer, or X, Y, Z in the 
case of a plane array of transducers), at least one 
shear wave propagation parameter is calculated, either 
at certain points of the observation field which are 
chosen by the user from the microcomputer 4, or 
throughout the observation field. 

The shear wave propagation parameter which is 
calculated in the course of the mapping step is chosen 
for example from the speed Cs of the shear waves, the 
shear modulus ^, the attenuation a of the shear waves, 
the shear elasticity jxl, shear viscosity /x2 . 

This calculation is performed via a conventional 
inversion process, an example of which is given 
hereinbelow in the case of a two-dimensional 
observation field (the same process would apply mutatis 
mutandis in the case of an observation field having one 
or three dimensions, respectively for a single 
transducer Tl or for a plane array of transducers) . 

In this example, the approximation will be made that 
the shear viscosity y.2 is zero and that the medium is 
isotropic . 

The wave equation giving the displacement vector D of 
each point of the medium 1 may be written: 

d 2 D — 
p— = A(tiJ)) (I), 

dr 

where p is the density of the medium 1, fi is the shear 
modulus (by assumption reduced to its real part the 
shear elasticity /xl since the shear viscosity /i2 is 
assumed to be zero) . 
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For the first component u of the vector D, that is to 
say for the displacement of the medium 1 in the 
direction X, we therefore have: 

5 p ^'if + -^ 1111 • 

After temporal Fourier transform and discretization, 
this equation can be written in the following matrix 
form, which can be written for each frequency of the 
10 spectrum of the shear wave: 

B=HM (III) , where: 

M is a vector of dimension (L+2) . (M+2) -4 , each 

15 component of which equals Mim, that is to say the local 
value of the shear modulus at each discretized point 
with coordinates (xi,y m ), where 1 and m are integers 
lying between 0 and respectively L+l and M+l, 
eliminating the pairs (l,m) equal to (0,0), (0,M+1), 

20 (L+1,0) and (L+l, M+l), L+2 and M+2 being the numbers of 
discretized points in the image of the medium 1 
respectively along the X and Y axes, 

B is a vector of dimension L.M, whose 
components equal -tij 2 .p.Ui m where w is the angular 

25 frequency of the low- frequency shear wave, p is the 
density of the medium, U im is the temporal Fourier 
transform of the displacement u at the point with 
coordinates (x^yj, 1 lying between 1 and L and m lying 
between 1 and M, 

30 - and H is a matrix of dimension L.M rows by 

(L+2) . (M+2) -4 columns, all of whose components are also 
known from the wave equation. 

By juxtaposing a sufficient number of equations (III) 
3 5 corresponding respectively to various available items 
of the frequency spectrum of the shear wave, a global 
matrix equation is obtained which can be solved by 
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matrix inversion, to obtain the vector M, that is to say the value of the shear 
modulus n at every point of the observation field. 


The local value of the speed of propagation Cs of the shear wave at 
5 each point can then be deduced therefrom, if desired, through the formula: 



10 

The mode of calculation would be the same if using the strains of the 
observed medium 1 rather than the displacements. 

Throughout the description and claims of this specification, the word 
15 "comprise' 5 and variations of that word, such as "comprising" and "comprises" 
are not intended to exclude other additives, steps or integers. 


• * • 
• 
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THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS 
1. A method of imaging using shear waves to observe a diffusing 
viscoelastic medium which contains particles reflecting ultrasonic 
compression waves, in which method an elastic shear wave is generated in the 
5 viscoelastic medium and the displacement of the viscoelastic medium 
subjected to said shear wave is observed by means of at least one ultrasonic 
compression wave, 

wherein the shear wave is generated by applying to the viscoelastic 
medium an excitation having the form of a low-frequency pulse which 
10 exhibits a central frequency f of between 20 and 5000 Hz, this low-frequency 
pulse exhibiting a duration of between l/2f and 20/f, 

wherein it comprises a propagation observation step in the course of 
which the propagation of the shear wave is observed simultaneously at a 
multitude of points in the observed medium, these points forming a 
1 5 substantially continuous observation field extending at least along a first axis, 
this shear wave propagation observation step consisting in: 

- emitting into the observed medium a succession of at least 10 shots of 
ultrasonic compression waves at a rate of between 100 and 100 000 shots per 
second, 

20 - detecting and recording in real time the echoes generated by the 

reflecting particles of the viscoelastic medium at each ultrasonic wave shot, 
these echoes corresponding to successive images of the observed medium, and 
wherein said method furthermore comprises a subsequent image 
processing step in the course of which the images thus obtained are processed 

25 at a later time at least by cross-correlation between successive images, so as to 
determine at each point of the observation field a motion parameter choserf 
between the displacement and the strain of the viscoelastic medium, in such a 
way as to thus obtain a succession of images showing the evolution of the 
motion parameter of the viscoelastic medium under the effect of the 

30 propagation of the shear wave. 
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2. The method as claimed in claim 1, in which the duration of the low- 
frequency pulse is between l/2f and 2f. 

3. The method as claimed in claim 1 or claim 2, in which the central 
frequency of the low-frequency pulse is between 30 and 1000 Hz. 

4. The method as claimed in any one of the preceding claims, in which the 
observed viscoelastic medium consists of a living body comprising at least one 
internal organ subjected to pulsatile motions, the low-frequency pulse which 
generates the shear wave being constituted by a pulsatile motion of said 
internal organ. 

5. The method as claimed in any one of claims 1 to 3, in which the 
observed viscoelastic medium is delimited by an outside surface and the low- 
frequency pulse is applied in the vicinity of this outside surface. 

6. The method as claimed in claim 5, in which the low-frequency pulse is 
applied by a means of excitation chosen from: 

- an acoustic wave generated by at least one acoustic transducer, and 

- a shock generated locally by physical contact in the vicinity of the 
outside surface of the viscoelastic medium. 

7. The method as claimed in claim 5 or claim 6, in which the ultrasonic 
compression wave shots are emitted and the echoes generated by the reflecting 
particles of the viscoelastic medium are detected by means of a bank of 
transducers which comprises at least one transducer and which is arranged in 
contact with the outside surface of the viscoelastic medium, the shear wave 
being applied to the viscoelastic medium by imposing a pulsatile displacement 
on said bank of transducers. 
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8. The method as claimed in claim 7, in which said motion parameter is 
the strain of the viscoelastic medium. 


9. The method as claimed in any one of the preceding claims in which, in 
5 the course of the observation of the propagation of the shear wave, between 

100 and 10 000 ultrasonic compression wave shots are emitted at a rate of 
between 1000 and 100 000 shots per second. 

10. The method as claimed in any one of the preceding claims, in which the 
1 0 observation field extends at least along a plane comprising on the one hand, 

the first axis and on the other hand, a second axis perpendicular to the first 
axis. 

11. The method as claimed in claim 10, in which, in the course of the 
15 propagation observation step, a bank of several acoustic transducers arranged 

at least along the second axis is used to emit the ultrasonic acoustic 
compression wave shots and detect the echoes generated by the reflecting 
particles of the viscoelastic medium, the echoes detected by each acoustic 
transducer being stored directly without prior processing in the course of the 
20 propagation observation step, and the image processing step comprising a 
preliminary substep of forming paths in the course of which an image of the 
viscoelastic medium corresponding to each ultrasonic compression wave shot 
is generated by combining at least some of the echoes received by the various 
transducers. 

25 

12. The method as claimed in any one of the preceding claims, in which the 
image processing step is followed by a viewing step in the course of which a 
film consisting of the succession of processed images is viewed under slow 
motion, each point of each image exhibiting an optical parameter which varies 

30 according to the value of the motion parameter assigned to this point. 
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13. The method as claimed in claim 12, in which the optical parameter is 
chosen from the gray level and the chromatic level. 


14. The method as claimed in any one of the preceding claims, in which the 
5 image processing step is followed by a mapping step in the course of which, 
on the basis of the evolution of the motion parameter over the course of time 
in the observation field, at least one propagation parameter is calculated for 
the shear wave at at least some points of the observation field. 

10 15. The method as claimed in claim 14, in which the shear wave 

• • • 

I propagation parameter which is calculated in the course of the mapping step is 

• • • 

chosen from the speed of the shear waves, the shear modulus, the attenuation 

« • • 

of the shear weaves, the shear elasticity, and the shear viscosity. 

■ • 

• • • • 

• • • 

15 16. An imaging device using shear waves to observe a diffusing 

, viscoelastic medium which contains particles reflecting ultrasonic 

• « • 

compression waves, this device comprising means of excitation for generating 
an elastic shear wave in the viscoelastic medium and means of acquisition for 

•• • 

observing, by means of at least one ultrasonic compression wave, the 
20 displacement of the viscoelastic medium subjected to said-shear wave, 

• • • • 

• wherein the means of excitation are adapted to apply to the viscoelastic 

medium an excitation having the form of a low-frequency pulse which 
exhibits a central frequency f of between 20 and 5000 Hz, this low-frequency 
pulse exhibiting a duration of between l/2f and 20/f, and 
25 wherein the means of acquisition are adapted to observe the 

propagation of the shear wave simultaneously at a multitude of points in the 
observed medium, these points forming a substantially continuous observation 
field extending at least along a first axis, said means of acquisition being 
adapted for: 


• • • 

• • • 

• • 

• ••• 

• • • 

• • • 

• • • 

• • • 

• • • 


» ■ 


• ••• 

• • 

• • • 
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- emitting into the observed medium a succession of at least 10 shots of 
ultrasonic compression waves at a rate of between 100 and 100 000 shots per 
second, 

- detecting and recording in real time the echoes generated by the 
5 reflecting particles of the viscoelastic medium at each ultrasonic wave shot, 

these echoes corresponding to successive images of the observed medium, and 
wherein said device furthermore comprises image processing means 
adapted to process at a later time the images obtained by the observation 
means, at least by cross-correlation between successive images, so as to 
10 determine at each point of the observation field a motion parameter chosen 
between the displacement and the strain of the viscoelastic medium, in such a 
way as to thus obtain a succession of images showing the evolution of the 
motion parameter of the viscoelastic medium under the effect of the 
propagation of the shear wave. 

15 

17. The device as claimed in claim 16, in which the means of observation 
comprise a bank of transducers which includes at least one transducer and 
which is adapted to be arranged in contact with an outside surface delimiting 
the viscoelastic medium, the means of excitation being adapted to impose a 

20 pulsatile displacement on said bank of transducers. 

18. The device as claimed in claim 17, in which said motion parameter is 
the strain of the viscoelastic medium. 


25 19. The device as claimed in any one of claims 17 and 18, in which the 
observation field extends at least along a plane comprising on the one hand, 
the first axis and on the other hand, a second axis perpendicular to the first 
axis, the bank of transducers comprising several transducers arranged at least 
along the second axis, control means being provided for selectively operating 

30 the device either in the mode of imaging by shear waves, or in a standard 
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echography mode making it possible to acquire between 10 and 100 images 
per second. 

20. A method of imaging using shear waves substantially as hereinbefore 
5 described with reference to the accompanying drawings. 

21. An imaging device using shear waves substantially as hereinbefore 
described with reference to the accompanying drawings. 

10 DATED: 23 September, 2002 
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